Admixture of divergent genomes facilitates hybridization across species in the family Brassicaceae.

Hybridization and polyploidization are pivotal to plant evolution. Genetic crosses between distantly related species are rare in nature due to reproductive barriers but how such hurdles can be overcome is largely unknown. Here we report the hybrid genome structure of xBrassicoraphanus, a synthetic allotetraploid of Brassica rapa and Raphanus sativus. We performed cytogenetic analysis and de novo genome assembly to examine chromosome behaviors and genome integrity in the hybrid. Transcriptome analysis was conducted to investigate expression of duplicated genes in conjunction with epigenome analysis to address whether genome admixture entails epigenetic reconfiguration. Allotetraploid xBrassicoraphanus retains both parental chromosomes without genome rearrangement. Meiotic synapsis formation and chromosome exchange are avoided between nonhomologous progenitor chromosomes. Reconfiguration of transcription network occurs, and less divergent cis-elements of duplicated genes are associated with convergent expression. Genome-wide DNA methylation asymmetry between progenitors is largely maintained but, notably, B. rapa-originated transposable elements are transcriptionally silenced in xBrassicoraphanus through gain of DNA methylation. Our results demonstrate that hybrid genome stabilization and transcription compatibility necessitate epigenome landscape adjustment and rewiring of cis-trans interactions. Overall, this study suggests that a certain extent of genome divergence facilitates hybridization across species, which may explain the great diversification and expansion of angiosperms during evolution.

Universal gene co-expression network reveals receptor-like protein genes involved in broad-spectrum resistance in pepper (Capsicum annuum L.).

Receptor-like proteins (RLPs) on plant cells have been implicated in immune responses and developmental processes. Although hundreds of RLP genes have been identified in plants, only a few RLPs have been functionally characterized in a limited number of plant species. Here, we identified RLPs in the pepper (Capsicum annuum) genome and performed comparative transcriptomics coupled with the analysis of conserved gene co-expression networks (GCNs) to reveal the role of core RLP regulators in pepper-pathogen interactions. A total of 102 RNA-seq datasets of pepper plants infected with four pathogens were used to construct CaRLP-targeted GCNs (CaRLP-GCNs). Resistance-responsive CaRLP-GCNs were merged to construct a universal GCN. Fourteen hub CaRLPs, tightly connected with defense-related gene clusters, were identified in eight modules. Based on the CaRLP-GCNs, we evaluated whether hub CaRLPs in the universal GCN are involved in the biotic stress response. Of the nine hub CaRLPs tested by virus-induced gene silencing, three genes (CaRLP264, CaRLP277, and CaRLP351) showed defense suppression with less hypersensitive response-like cell death in race-specific and non-host resistance response to viruses and bacteria, respectively, and consistently enhanced susceptibility to Ralstonia solanacearum and/or Phytophthora capsici. These data suggest that key CaRLPs are involved in the defense response to multiple biotic stresses and can be used to engineer a plant with broad-spectrum resistance. Together, our data show that generating a universal GCN using comprehensive transcriptome datasets can provide important clues to uncover genes involved in various biological processes.

Prometheus, an omics portal for interkingdom comparative genomic analyses.

Functional analyses of genes are crucial for unveiling biological responses, genetic engineering, and developing new medicines. However, functional analyses have largely been restricted to model organisms, representing a major hurdle for functional studies and industrial applications. To resolve this, comparative genome analyses can be used to provide clues to gene functions as well as their evolutionary history. To this end, we present Prometheus, a web-based omics portal that contains more than 17,215 sequences from prokaryotic and eukaryotic genomes. This portal supports interkingdom comparative analyses via a domain architecture-based gene identification system and Gene Search, and users can easily and rapidly identify single or entire gene sets in specific pathways. Bioinformatics tools for further analyses are provided in Prometheus or through Bio-Express, a cloud-based bioinformatics analysis platform. Prometheus is a new paradigm for comparative analyses of large amounts of genomic information.

TGFam-Finder: a novel solution for target-gene family annotation in plants.

Whole-genome annotation error that omits essential protein-coding genes hinders further research. We developed Target Gene Family Finder (TGFam-Finder), an alternative tool for the structural annotation of protein-coding genes containing target domain(s) of interest in plant genomes. TGFam-Finder took considerably reduced annotation run-time and improved accuracy compared to conventional annotation tools. Large-scale re-annotation of 50 plant genomes identified an average of 150, 166 and 86 additional far-red-impaired response 1, nucleotide-binding and leucine-rich-repeat, and cytochrome P450 genes, respectively, that were missed in previous annotations. We detected significantly higher number of translated genes in the new annotations using mass spectrometry data from seven plant species compared to previous annotations. TGFam-Finder along with the new gene models can provide an optimized platform for comprehensive functional, comparative, and evolutionary studies in plants.

Comprehensive analysis of Translationally Controlled Tumor Protein (TCTP) provides insights for lineage-specific evolution and functional divergence.

BACKGROUND: Translationally controlled tumor protein (TCTP) is a conserved, multifunctional protein involved in numerous cellular processes in eukaryotes. Although the functions of TCTP have been investigated sporadically in animals, invertebrates, and plants, few lineage-specific activities of this molecule, have been reported. An exception is in Arabidopsis thaliana, in which TCTP (AtTCTP1) functions in stomatal closuer by regulating microtubule stability. Further, although the development of next-generation sequencing technologies has facilitated the analysis of many eukaryotic genomes in public databases, inter-kingdom comparative analyses using available genome information are comparatively scarce. METHODOLOGY: To carry out inter-kingdom comparative analysis of TCTP, TCTP genes were identified from 377 species. Then phylogenetic analysis, prediction of protein structure, molecular docking simulation and molecular dynamics analysis were performed to investigate the evolution of TCTP genes and their binding proteins. RESULTS: A total of 533 TCTP genes were identified from 377 eukaryotic species, including protozoa, fungi, invertebrates, vertebrates, and plants. Phylogenetic and secondary structure analyses reveal lineage-specific evolution of TCTP, and inter-kingdom comparisons highlight the lineage-specific emergence of, or changes in, secondary structure elements in TCTP proteins from different kingdoms. Furthermore, secondary structure comparisons between TCTP proteins within each kingdom, combined with measurements of the degree of sequence conservation, suggest that TCTP genes have evolved to conserve protein secondary structures in a lineage-specific manner. Additional tertiary structure analysis of TCTP-binding proteins and their interacting partners and docking simulations between these proteins further imply that TCTP gene variation may influence the tertiary structures of TCTP-binding proteins in a lineage-specific manner. CONCLUSIONS: Our analysis suggests that TCTP has undergone lineage-specific evolution and that structural changes in TCTP proteins may correlate with the tertiary structure of TCTP-binding proteins and their binding partners in a lineage-specific manner.

Transcriptome profiling of abiotic responses to heat, cold, salt, and osmotic stress of Capsicum annuum L.

Peppers (Capsicum annuum L.), belonging to the Solanaceae family, are one of the most economically important crops globally. Like other crops, peppers are threatened by diverse environmental conditions due to different pathogens and abiotic stresses. High-quality reference genomes with massive datasets of transcriptomes from various conditions can provide clues to preferred agronomic traits for breeding. However, few global gene expression profiling datasets have been published to examine the environmental stress-resistant mechanisms in peppers. In this study, we report the RNA-seq analyses of peppers treated with heat, cold, salinity, and osmotic stress at six different time points. RNA-seq libraries from 78 RNA samples containing three biological replicates per time point for each of the abiotic stresses and a mock control were constructed. A total of 204.68 Gb of transcriptome data were verified by differentially expressed genes and gene ontology enrichment analysis. Analyses of the transcriptome data in this study will provide useful information for basic studies of various stimuli to facilitate the development of stress-resistant pepper cultivars.

Draft genome sequences of two oriental melons, Cucumis melo L. var. makuwa.

Oriental melon (Cucumis melo L. var. makuwa) is one of the most important cultivated cucurbits, and is grown widely in Northeast Asian countries. With increasing interest in its biological properties and economic importance, oriental melon has become an attractive model crop for studying various horticultural traits. A previous genome sequence of the melon was constructed from a homozygous double-haploid line. Thus, individual reference genomes are required to perform functional studies and further breeding applications. Here, we report draft genome sequences of two oriental melons, Chang Bougi and SW3. The assembled 344 Mb genome of Chang Bougi was obtained with scaffold N50 1.0 Mb, and 36,235 genes were annotated. The 354 Mb genome of SW3 was assembled with scaffold N50 1.6 Mb, and has 38,173 genes. These newly constructed genomes will enable studies of fruit development, disease resistance, and breeding applications in the oriental melon.

Genome Sequence of Striga asiatica Provides Insight into the Evolution of Plant Parasitism.

Parasitic plants in the genus Striga, commonly known as witchweeds, cause major crop losses in sub-Saharan Africa and pose a threat to agriculture worldwide. An understanding of Striga parasite biology, which could lead to agricultural solutions, has been hampered by the lack of genome information. Here, we report the draft genome sequence of Striga asiatica with 34,577 predicted protein-coding genes, which reflects gene family contractions and expansions that are consistent with a three-phase model of parasitic plant genome evolution. Striga seeds germinate in response to host-derived strigolactones (SLs) and then develop a specialized penetration structure, the haustorium, to invade the host root. A family of SL receptors has undergone a striking expansion, suggesting a molecular basis for the evolution of broad host range among Striga spp. We found that genes involved in lateral root development in non-parasitic model species are coordinately induced during haustorium development in Striga, suggesting a pathway that was partly co-opted during the evolution of the haustorium. In addition, we found evidence for horizontal transfer of host genes as well as retrotransposons, indicating gene flow to S. asiatica from hosts. Our results provide valuable insights into the evolution of parasitism and a key resource for the future development of Striga control strategies.

Identification of transcriptome-wide, nut weight-associated SNPs in Castanea crenata.

Nut weight is one of the most important traits that can affect a chestnut grower's returns. Due to the long juvenile phase of chestnut trees, the selection of desired characteristics at early developmental stages represents a major challenge for chestnut breeding. In this study, we identified single nucleotide polymorphisms (SNPs) in transcriptomic regions, which were significantly associated with nut weight in chestnuts (Castanea crenata), using a genome-wide association study (GWAS). RNA-sequencing (RNA-seq) data were generated from large and small nut-bearing trees, using an Illumina HiSeq. 2000 system, and 3,271,142 SNPs were identified. A total of 21 putative SNPs were significantly associated with chestnut weight (false discovery rate [FDR] < 10-5), based on further analyses. We also applied five machine learning (ML) algorithms, support vector machine (SVM), C5.0, k-nearest neighbour (k-NN), partial least squares (PLS), and random forest (RF), using the 21 SNPs to predict the nut weights of a second population. The average accuracy of the ML algorithms for the prediction of chestnut weights was greater than 68%. Taken together, we suggest that these SNPs have the potential to be used during marker-assisted selection to facilitate the breeding of large chestnut-bearing varieties.

New reference genome sequences of hot pepper reveal the massive evolution of plant disease-resistance genes by retroduplication.

BACKGROUND: Transposable elements are major evolutionary forces which can cause new genome structure and species diversification. The role of transposable elements in the expansion of nucleotide-binding and leucine-rich-repeat proteins (NLRs), the major disease-resistance gene families, has been unexplored in plants. RESULTS: We report two high-quality de novo genomes (Capsicum baccatum and C. chinense) and an improved reference genome (C. annuum) for peppers. Dynamic genome rearrangements involving translocations among chromosomes 3, 5, and 9 were detected in comparison between C. baccatum and the two other peppers. The amplification of athila LTR-retrotransposons, members of the gypsy superfamily, led to genome expansion in C. baccatum. In-depth genome-wide comparison of genes and repeats unveiled that the copy numbers of NLRs were greatly increased by LTR-retrotransposon-mediated retroduplication. Moreover, retroduplicated NLRs are abundant across the angiosperms and, in most cases, are lineage-specific. CONCLUSIONS: Our study reveals that retroduplication has played key roles for the massive emergence of NLR genes including functional disease-resistance genes in pepper plants.

Transcriptome analysis of the oriental melon (Cucumis melo L. var. makuwa) during fruit development.

BACKGROUND: The oriental melon (Cucumis melo L. var. makuwa) is one of the most important cultivated cucurbits grown widely in Korea, Japan, and northern China. It is cultivated because its fruit has a sweet aromatic flavor and is rich in soluble sugars, organic acids, minerals, and vitamins. In order to elucidate the genetic and molecular basis of the developmental changes that determine size, color, and sugar contents of the fruit, we performed de novo transcriptome sequencing to analyze the genes expressed during fruit development. RESULTS: We identified a total of 47,666 of representative loci from 100,875 transcripts and functionally annotated 33,963 of the loci based on orthologs in Arabidopsis thaliana. Among those loci, we identified 5,173 differentially expressed genes, which were classified into 14 clusters base on the modulation of their expression patterns. The expression patterns suggested that the differentially expressed genes were related to fruit development and maturation through diverse metabolic pathways. Analyses based on gene set enrichment and the pathways described in the Kyoto Encyclopedia of Genes and Genomes suggested that the expression of genes involved in starch and sucrose metabolism and carotenoid biosynthesis were regulated dynamically during fruit development and subsequent maturation. CONCLUSION: Our results provide the gene expression patterns related to different stages of fruit development and maturation in the oriental melon. The expression patterns give clues about important regulatory mechanisms, especially those involving starch, sugar, and carotenoid biosynthesis, in the development of the oriental melon fruit.

Genome analysis of Hibiscus syriacus provides insights of polyploidization and indeterminate flowering in woody plants.

Hibiscus syriacus (L.) (rose of Sharon) is one of the most widespread garden shrubs in the world. We report a draft of the H. syriacus genome comprised of a 1.75 Gb assembly that covers 92% of the genome with only 1.7% (33 Mb) gap sequences. Predicted gene modeling detected 87,603 genes, mostly supported by deep RNA sequencing data. To define gene family distribution among relatives of H. syriacus, orthologous gene sets containing 164,660 genes in 21,472 clusters were identified by OrthoMCL analysis of five plant species, including H. syriacus, Arabidopsis thaliana, Gossypium raimondii, Theobroma cacao and Amborella trichopoda. We inferred their evolutionary relationships based on divergence times among Malvaceae plant genes and found that gene families involved in flowering regulation and disease resistance were more highly divergent and expanded in H. syriacus than in its close relatives, G. raimondii (DD) and T. cacao. Clustered gene families and gene collinearity analysis revealed that two recent rounds of whole-genome duplication were followed by diploidization of the H. syriacus genome after speciation. Copy number variation and phylogenetic divergence indicates that WGDs and subsequent diploidization led to unequal duplication and deletion of flowering-related genes in H. syriacus and may affect its unique floral morphology.